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Magnetization and energy gap of Kondo insulator YbBi2 are calculated theoretically based on 
the previously proposed tight-binding model composed of Yb 5de and 4f T$ orbitals. It is found 
that magnetization curves are almost isotropic, naturally expected from the cubic symmetry, 
but that the gap-closing field has an anisotropy: the gap closes faster for the field in (100) 
direction than in (110) and (111) directions, in accord with the experiments. This is qualitatively 
understood by considering the maximal eigenvalues of the total angular momentum operators 
projected on each direction of the magnetic field. But the numerical calculation based on the 
band model yields better agreement with the experiment. 
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1. Introduction 

YbBi2 X is the most typical Kondo insulator, since it 
has a cubic crystal structure (Yb and B12 forms NaCl 
structure), and many experiments: transport, 2 calorimet- 
ric, 3 photoemission, 4 ' 5 optical, 6,7 neutron scattering, 8-10 
etc. are already reported. We believe that most of them 
can be understood in terms of the band insulator model 
with strong Coulomb repulsion between f electrons. In 
fact, based on the LDA+U band calculation, it was found 
that the conduction bands near the Fermi level are com- 
posed of the 5de orbitals on Yb atoms, and the simple 
tight-binding band reproduces the LDA+U calculation 
rather well. 11 Some of the conduction bands are dou- 
bly degenerate (four- fold if spin degeneracy is included) . 
A gap opens by the hybridization of the four-fold de- 
generate T% states with the above-mentioned conduction 
bands with degeneracy. It was pointed out that the gap 
can not open if the crystal-field groud state of 4f were 
doublet. Therefore, proper consideration of the degener- 
acy is a key to understanding the opening of a gap in the 
Kondo insulators. 11 

Based on this band model, one of the authors 12 calcu- 
lated frequency- and temperature-dependence of the op- 
tical conductivity and obtained semi-quantitative agree- 
ment with the experiment. 7 He succeeded in reproducing 
the mid-infrared peak and the prominent temperature 
dependence of the low frequency part, which is due to 
the many-body effect. 

In contrast to YbBi2, it is not clear whether other 
well-known or classical "Kondo insulators", SmS 13 and 
TmSe, 14 belong to the same category as YbBi2, since 
Sm and Tm has more complicated 4f states than Yb. In 
addition, TmSe has odd number of electrons per unit cell, 
so that it can not be a band insulator if 4f electrons are 
mixed valent. Other Kondo insulators, which may belong 
to the same category as YbBi2, will be Ce3Bi4Pt3 15 and 
Ce-skutterudites, e.g. CeFe4Asi2- 16 Nevertheless, YbBi2 



can be a prototype of the Kondo insulator from the point 
of view of the simpleness of the band structure, sufficient 
number of experiments and the existence of the compact 
theoretical model. 

Despite the above-mentioned success, YbBi2 still has 
many features to be understood correctly and quantita- 
tively. One of them is the magnetic properties. Sugiyama, 
et al. 2 first found that the activation energy of the re- 
sistivity of this insulator vanishes at aroud 50 Tesla 
when magnetic field is applied. This is called as field- 
induced insulator-to-metal transition and the critical 
field is denoted as B c hereafter. Iga, et al. 17 found that 
the magnetization curve for each direction does not show 
an anisotropy, but the critical field B c does. Namely, 
B c (110) and S c (lll) is larger than B c (100) by 10 per- 
cent (B c (100) =47T, B c (110) ~ B c (lll) =53T). The 
transitions are first order since hysterisis is observed. 

Using the simple periodic Anderson model, one of the 
present authors showed that the overall feature of the 
magnetization curve of YbB 12 can be explained and the 
transition can be first order if the Coulomb repulsion is 
taken into account. 18, 19 In addition, he found that the 
renormalization factor (which renormalizes the gap) due 
to strong correlation does not largely change by the mag- 
netic field in the insulating phase, so that the renormal- 
ized up- and down-spin bands are shifted by the field as 
if they are rigid bands and the renormalized band gap is 
closed also as if in the case of the rigid bands. 18 But the 
anisotropic feature of the f-electrons under cubic crys- 
talline field was not taken into account. 

In the present study, we investigate theoretically the 
origin and mechanism of the anisotropy of the field- 
induced insulator-to-metal transition. In §2, we briefly 
explain the previously introduced tight-binding band 
model. Calculations of magnetization and energy gap un- 
der magnetic field will be presented in §3. The last section 
is devoted to the discussions and conclusions. 
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2. Tight-Binding Energy Band 

The energy dispersions of the conduction bands of 
YbBi2 are well approximated by the tight-binding model 
of 5de orbitals of Yb on the fee lattice and expressed by 
the following simple expressions: 11 

K P = E de + 3(dda) cos(fc Q /2) cos(fc^/2), (1) 

where (ddc) is a Slater-Koster integral, 20 (a, (3) denotes 
(x, y), (y, z) or (z, x). The lowest band along r-K-X(llO) 
and r-X(lOO) are doubly degenerate (four-fold if the spin 
degeneracy is included). Introducing the hybridization 
between these d bands and the 4f j = 7/2 T$ states, we 
obtained the energy bands with a gap. 

In eq.(l), we locate the energy level of 5de orbitals at 
E d£ =1.0 Ry, and set (ddcr)= 0.06 Ry. This value of (ddcr) 
is the effective one due to hopping between near-neighbor 
Yb sites through B 12 clusters. 

The hybridization between these 5d bands and the 4f 
states is described by the effective (dfcx) integrals 21 be- 
tween the nearest-neighbor Yb sites. We retain only the 
nearest neighbor (dfer) bonds as the simplest model. We 
locate the 4f T$ states at Er s = 0.88 Ry and choose 
(dfo-)=0.015 Ry. We also include (df?r) = -0.0075 Ry 
and (ffcr) = -0.003 Ry to adjust to the LDA+U cal- 
culation. Furthermore, the filled bands below the gap 
are shifted down by AE = —0.011 Ry relative to the 
bands above the gap. This treatment is in accord with 
the spirit of the LDA+U treatment. These parameters 
are the same as those in our previous paper 12 to fit the 
optical conductivity. 7 

Using these parameters, we obtain the dispersion 
curves which reproduce the LDA+U bands around the 
gap rather well. 11 They have an indirect gap of about 
0.0069 Ry (= 1089 K) between X and L points and the 
direct gap of 0.018 Ry (= 0.245 eV). The former is about 
six times larger than the best estimated value of the in- 
direct gap in the optical experiment, E g =174 K. 7 How- 
ever, in the present realistic band model, we can not use 
this value. If we use a band model with the correct value 
of the indirect gap, it is difficult to reproduce the gap in 
the density of states, since we must use a very fine mesh 
to calculate the three-dimensional fc-summation over the 
Brillouin zone. 

Note that the band calculation can treat the f-state 
only by the j-j coupling scheme. One hole in j = 7/2 T$ 
state corresponds to the T 8 state with the total angu- 
lar momentum J = 7/2 of Yb 3+ (the f-electron number 
nf — 13). In the present case, the number of filled 4f 
r 8 electrons is about 3.1 and 0.9 hole exists in the bands 
above the gap, whereas j — 5/2, j = 7/2 Tq and T-? states 
are far below Tg and are implicitly treated as completely 
filled. Thus the present model approximately reproduces 
the experimental valence of Yb +2 9 . 

3. Calculation of Magnetization and Energy 
Gap under Magnetic Field 

We apply the magnetic field B = fi H along (100), 
(110) and (111) directions. The Hamiltonian reads 
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Fig. 1. The tight-binding band for YbBi2 under B=500T in (100) 
direction. The X points in the left and right denote the equivalent 
(110) and (100) points (in unit of 2n/a), respectively. 



U = ^2 t tf c L c 30+Yl3c^BS l -H + Y^9.JHBji-H, (2) 

a, i i 

i,j 

where denotes the hopping matrix elements between 
the ten orbitals {a, j3 = xy \,yz \,zx \,xy l,yz J., zx J. 
, Tg 1 ^, Tg 2 ^} on i and j sites and is expressed by the 
Slater-Koster parameter given in the last section. 12 /j,b 
is the Bohr magneton, g c = 2 and gf = 8/7 are the g- 
factors of conduction and f electrons. Sj and j i are the 
spin and the total angular momentum operators for the 
conduction and f electron at i site, respectively. We have 
omitted the effect of magnetic field on the orbital motion 
except that through j. The effect of orbital motion can 
be included by attaching the phase factor to the hopping 
matrix elements, but it will be smaller than those con- 
sidered here through j, since the latter will be enhanced 
by the strong correlation. 

In the present model with the larger gap than the ex- 
periment, the gap closes at the large field of about 500 
T. We show the band structure at this field in (100) di- 
rection in Fig.l. Four bands lie below the gap and six 
bands above the gap and all the degeneracies under zero 
field are lifted. 

Diagonalizing the Hamiltonian, the Bloch function is 
obtained as 

ipkn{r) = --L «kn, Q e tk - R '^(r - R^, (3) 

where 4> a (r — Ri) is the local orbital a at the site Ri, 
n denotes the diagonalizcd band index, Ukn,a the coef- 
ficient and N the total number of unit cell. Using this, 
the magnetization M is calculated by 

M = ^^/(£kn)[0cMB<Vw|S'|V'k„} 
k.n 

+ ff/MB(^kn|J|^kn)], (4) 

where S — J2 i Sj, J = J2iJii f(E) ^ s the Fermi function 
and i?kn the n-th band energy. Note that the applica- 
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tion of the magnetic field lowers the symmetry, so that 
the sum over 1/48 of the first Brillouine zone is not al- 
lowed. In the actual calculation, we summed over whole 
of the first and the second Brillouine zones because the 
programming becomes much easier: we can simply make 
summation over a cube with twice the volume of the first 
Brillouine zone of fee lattice. We divide this cube into at 
most 400 3 mesh points. Temperature is set equal to zero 
throughout the paper. 

Also for simplicity, we have summed only over the low- 
est four filled bands in the calculation of the magnetiza- 
tion, so that the calculation makes sense only for B < B c 
(critical field) in each field direction. For B > B c , one has 
to determine the chemical potential for each value and 
direction of the field and sum up over more than four 
bands. We did not do it since we are interested only in 
the anisotropy problem in the present paper. 

M/hb = \M\/fj,B is plotted in Fig. 2 for the three di- 
rections of the field. In accord with the general theorem 
of the group theory that the linear response tensor must 
be proportional to the unit tensor (isotropic), the present 
result shows that the magnetization is isotropic up to 200 
T. (Indirect) energy gap vs. magnetic field is plotted in 
Fig. 3. Since we take large value of the initial gap, the 
critical value of the field in each direction is also large: 
B c (100) = 464T, B c (110) = 501T, B c (lll) = 506T. 
These values will become smaller if we use the correct 
value of the band gap. For B < 300T, E g (lll) becomes 
slightly smaller than i? s (110), but the difference is very 
small and it seems out of range of the precision of the 
present calculation. It seems reasonable to consider that 
E g (100) < E g (U0)<E g (Ul) for all the field, as will bc 
discussed in the next section. 




B[T] 

Fig. 2. Magnetization curves under the magnetic field in (100 
(110) and (111) directions. 



4. Discussions and Conclusions 

In order to understand an origin of the anisotropy < 
the critical field, we define the angular momentum ope: 
ators projected onto each direction: j 100 = j ■ u wo = j, 
3 no = 3 ■ "no = (jx + j v )/V2, j ul = j ■ Mm = 
(jx +j y + .7z)/v / 3, where Uioo, "no and u m are the un 
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Fig. 3. Closing of the energy gap by the magnetic field in (100), 
(110) and (111) directions. Difference between (110) and (111) 
curves is almost invisible. 



vectors in each direction. Making the matrix representa- 
tion of these operators in the r§ space and diagonalizc 
the 4x4 matrices, we obtain the maximal eigenvalue for 
each direction as follows: 1.833 in (100), 1.635 in (110) 
and 1.500 in (111). From this result, one can understand 
that the f-states above and below the gap is shifted by the 
field faster in this order. (See Fig. 4.) There is, however, 
a difference between the above explanation and the ex- 
periment. Namely, the critical fields in (110) and (111) 
directions are almost the same, whereas i? c (100) is 10 
per cent smaller in the experiment. Our numerical cal- 
culation based on the realistic band model presented in 
Fig. 3. successfully reproduces this experimental feature. 
It may be because the percentage of each r 8 component 
is different at each k point, so that the above-mentioned 
explanation using the maximal eigenvalues of projected 
angular momentum operators holds only qualitatively. 
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Fig. 4. Schematic figure to explain how the gap is closed when 
the magnetic field is applied in each direction. 
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The other problem is that the anisotropy of the gap- 
closing is not clearly visible in Fig. 3 in the low field region 
where the magnetization is isotropic (B < 200 T). This 
may be partly because of insufficiency of the number of 
the fc-mesh points. We anticipate, however, that the clear 
anisotropy observed in the high field range in Fig. 3 is a 
real one and may survive also in the low field region in 
accord with the experiment. 

One more problem in our theory is that the calcu- 
lated magnetization is too small. Value of the magne- 
tization right before the metal-to-insulator transition is 
about 0.3 ~ 0.4 /Us/Yb in the experiment, 17 but the cal- 
culation shows only about 0.1 ^e/Yb even at 500 T, and 
M/fiB ~ 0.01 at H = 50 T. Namely, about 30-40 times 
enhancement is necessary. It is clear that the slope of the 
magnetization curve will be enhanced by a many-body 
effect, 18 but a quantitative calculation using the present 
model with full anisotropy is not easy. Nevertheless, it 
was already clarified that the renormarization factor does 
not strongly depend on the magnetic field until the gap 
is closed by the field at least for the simple periodic An- 
derson model 18 as mentioned in the Introduction. Since 
the magnetization must be isotropic in the linear region, 
the enhancement factor must also be isotropic. Thus, the 
many-body effect will give only a quantitative effect on 
the slopes of the magnetization curves and may not affect 
the anisotropy of the critical field. 

In conclusion, using the tight-binding model composed 
of Yb 5de and 4f r 8 orbitals, which reproduces the 
LDA+U band structure near the energy gap, we have cal- 
culated the magnetization curves and the collapse of the 
energy gap of the Kondo insulator YbBi2 when the mag- 
netic field is applied in (100), (110) and (111) directions. 
Since the initial gap was taken several times larger than 
the experimental one due to a technical reason in the 
numerical calculation, the critical values of the fields at 
which the gap vanishes are too large. We, however, have 
successfully explained the order of the direction of the 
critical field: B c (100) < S c (100) w B c (lll). Although 
many-body effect is not included, we have disscussed that 
the results may not change even if the strong correlation 
is taken into account. Of course, more perfect calculation 
should be performed including correlation and a small 
initial gap, but we believe that the present study can be 
an important step to understand the properties of the 
Kondo insulators from the microscopic point of view, es- 
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pecially, from the model based on the LDA(+U) band 
calculation plus correlation. 
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